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Abstract 

The breathing-mode isoscalar giant monopole resonance (GMR) is investigated us- 
ing the generator coordinate method within the relativistic mean-field (RMF) the- 
ory. Employing the Lagrangian models of the nonlinear-<7 model (NLcr), the scalar- 
vector interaction model (SVI) and the a-uj coupling model (SIGO), we show that 
each Lagrangian model exhibits a distinctly different GMR response. Consequently, 
Lagrangian models yield a different value of the GMR energy for a given value of the 
nuclear matter incompressibility K^. It is shown that this effect arises largely from 
a different value of the surface incompressibility K sur j inherent to each Lagrangian 
model, thus giving rise to the ratio K sur f jK^ which depends upon the Lagrangian 
model used. This is attributed to a difference in the density dependence of the me- 
son masses and hence to the density dependence of the nuclear interaction amongst 
various Lagrangian models. The sensitivity of the GMR energy to the Lagrangian 
model used and thus emergence of a multitude of GMR energies for a given value 
of -fToo renders the method of extracting K^ on the basis of interpolation amongst 
forces as inappropriate. As a remedy, the need to 'calibrate' the density dependence 
of the nuclear interaction in the RMF theory is proposed. 
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1 Introduction 



The compressibility of nuclear matter has been a matter of exploration and 
discussion for a long time. The incompressibility or the compression modulus 
is a fundamental property of the nuclear matter. It constitutes a cardinal point 
on the equation of state (EOS) of the nuclear matter. Defined as the second 
derivative of the EOS of the nuclear matter at the saturation density, it is 
also important for astrophysical phenomena such as supernova explosion and 
structure of neutron stars [Tf2"f5] . 

The breathing-mode isoscalar giant monopole resonance (GMR) is a principal 
source of information on the (in) compressibility of nuclei and nuclear matter 
[4llo1l6] . Nuclei undergo radial density oscillations about the equilibrium point 
in this mode. It is a small-amplitude collective motion wherein the nucleus 
participates as a whole. The isoscalar GMR is well-established experimentally 
and has been studied in a large number of nuclei from heavy to light mass re- 
gions in various laboratories around the world [7.8.9. 10J. Recent measurements 
|ll|fT2] have brought in improved precision in providing better information on 
the breathing-mode GMR. 

Experimental determination of the incompressibility of nuclear matter has 
entailed using a liquid-drop model type expansion of the incompressibility Ka 
of a nucleus in terms of a bulk and finite-size components [7|TU] . Experimental 
data on the breathing-mode isoscalar GMR on a large number of nuclei was 
used (7J. Attempts were made to extract the nuclear matter incompressibility 
Koo using the data on Sn and Sm isotopes obtained with considerable precision 
[TO] . However, the correlation between the Coulomb term and hinders 
extraction of the incompressibility of the nuclear matter [13J. 

With the difficulty of extracting directly using the experimental data, 
theoretical approach based upon interpolation amongst Skyrme and Gogny 
type interactions as proposed by Blaizot [4"][T4"] has been extensively employed. 
It has required reproducing experimental breathing-mode GMR energies of 
nuclei using an interaction with an appropriate value of in self-consistent 
Hartree-Fock and random-phase approximation (RPA) calculations [2]. In 
this scheme, the GMR energy of 208 Pb has been shown to be reproduced with 
forces with ~ 220 MeV. On the other hand, the RPA calculations with 
force(s) which reproduce the GMR energy of 208 Pb give much larger value for 
90 Zr, another key nucleus, as compared to the experimental data [14] . The 
GMR energy of 90 Zr [15] has since been overestimated by ~ 1 — 2 MeV in 
nonrelativistic [16] and relativistic [17] RPA approach. In the same vain, the 
GMR energies of Sn and Sm isotopes could not be described within the RPA 
sum-rule approach [T8"|[T9"] using the force SkM* with ~ 216 MeV, which 
described 208 Pb successfully though. It has been pointed out recently that 
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the GMR energies of Sn isotopes are overestimated in a self-consistent RPA 
approach by 1-2 MeV by the forces which reproduce the GMR energy of 208 Pb 
appropriately [2"D~|I2"T] . 

Using the RMF approach and relativistic RPA, a similar behaviour has been 
observed [22]. The success of the RPA approach and extraction of using 
the interpolation method thus hinges very strongly on the nucleus 208 Pb which 
masks attendant problems with other nuclei. 

With the advent of the relativistic mean-field (RMF) theory [2"3l24ll25|l2"o r ] . 
properties of nuclei and nuclear matter can be studied within this framework 
which has become successful in describing ground-state properties of nuclei 
along and far from the stability line [2711281129] . The inherent advantage of 
the RMF theory lies in the Dirac-Lorentz structure of nucleons, which allows 
a built-in spin-orbit interaction. Consequently, in properties related to shell 
structure, especially that of the anomalous behaviour of the isotope shifts of 
Pb nuclei, the RMF theory exhibits special advantages over the nonrelativistic 
Skyrme theory [30J. 

The breathing-mode isoscalar GMR has been studied within the RMF theory 
earlier [31]. Again, the theoretical or microscopic method of extracting Koq 
entails an interpolation amongst various forces with different values of 
with a view to being able to reproduce GMR energies of nuclei such as 208 Pb, 
Sn isotopes and 90 Zr among others within a microscopic approach [32] • The 
Lagrangian model with nonlinear scalar self-couplings of a meson has mostly 
been used with this approach. It has been pointed out in a schematic approach 
that forces with ~ 250 — 270 in the RMF model with the nonlinear scalar 
self-couplings (of the form <7 3 + cr 4 ) would be compatible with the experimental 
GMR energy of 208 Pb [32J. However, it has been shown recently that the 
GMR data on Sn isotopes and 90 Zr were overestimated by such a force |21j . 
Relativistic RPA calculations have also been performed with Lagrangians with 
density-dependent meson couplings [33], which focus only on the key nucleus 
of 208 Pb. 

The interpolation scheme has been the cornerstone of the theoretical method 
to extract the incompressibility of nuclear matter microscopically. Using the 
various Lagrangian models developed, it will be shown in this work that this 
approach runs into problems due to a sensitivity shown by GMR energies to 
the Lagrangian model employed. A difference in the density dependence of 
the meson masses amongst the various Lagrangian models and an ensuing 
difference in the response of the surface of nuclei to compression is ascribed 
for this behaviour. 

The subject matter of the paper is organized as follows: Definitions of the basic 
ingradients of the GMR and the incompressibility of the nuclear matter are 
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given in Section [2j In Section [3j the formalism of the RMF theory is provided. 
We discuss the Lagrangian models used in this work in Section HI Details of 
the generator coordinate method employed to explore and calculate the GMR 
energies of nuclei are discussed in Section It will be shown in Section E] that 
the GMR energies exhibit a sensitivity to the Lagrangian model considered. 
We shall show that the surface incompressibility plays an important role for 
the GMR energies in the RMF theory. The last section summarizes the main 
results and an outlook for future investigations is presented. 



2 The breathing-mode giant monopole resonance and the com- 
pressibility of nuclear matter 

The isoscalar (T=0) giant monopole resonance (GMR) has been well estab- 
lished in the last few decades [6]. It is a breathing-mode of a nucleus wherein 
the nucleus undergoes a density oscillation about its equilibrium value. The 
breathing mode GMR energy or its frequency is related to the incompressibil- 
ity Ka of a nucleus with mass number A by 



where m is the nucleon mass and (r 2 ) is the mean-square mass radius of the 
nucleus in its ground state. 

Analogous to the expansion of the total energy of a nucleus in terms of a 
liquid drop formula, one can write the expansion of the incompressibility of 
the nucleus in terms of contributions from the volume, surface, asymmetry, 
Coulomb and curvature terms, respectively, as [1]: 



where N and Z are the neutron and proton numbers, respectively. The K^s 
are 'incompressibilities' (coefficients) representing various terms and Cj's are 
the terms of the expansion in various powers of A. The volume term is 
defined as the curvature of the EOS of the infinite nuclear matter: 




(1) 



Ka = J2 KiCi 



+ K surf A-V* + K asym + K ( 




K, 




2 E/A 
dp 2 



(3) 
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where the second-derivative of the EOS, i.e., of e(p) is evaluated at the satu- 
ration density po- Thus, represents the most significant (reference) point 
of the high-density behaviour of the infinite nuclear matter. 

The volume and surface terms provide the most significant contributions to 
Ka, whereas the asymmetry and the Coulomb terms are much smaller in 
comparison. The curvature term is even smaller than the latter two. The 
volume term being the largest contributor to Ka influences the breathing- 
mode GMR energies of nuclei directly. 

The Coulomb term Kc ou i has been derived as [I] 



where e'" = d 3 (E/A)/dp 3 is the third-derivative (skewness) of the EOS. As 
readily seen from Eq. (j3J), Kcoui is correlated to K^. In Section [6] we shall see 
that the magnitude of the coefficient K Cou i is about two orders of magnitude 
smaller than the major contributors such as and K sur f. Due to this cor- 
relation superimposed on to the meagre quantity of the Coulomb term, it is 
difficult to extract from breathing- mode energies [13] . In the present work, 
we shall circumvent this correlation and attempt to disentangle the various 
terms. 



3 The Relativistic Mean-Field Theory - Formalism 

The starting point of the RMF theory is the basic Lagrangian (the linear 
Walecka model) that describes nucleons as Dirac spinors interacting with the 
meson fields [23] : 



where m is the bare nucleon mass and ip is its Dirac spinor. Nucleons in- 
teract with the a, u>, and p mesons. Here, g a , g u , and g p are the respective 
coupling constants of the interaction. The photonic field is produced by the 
electromagnetic vector A 11 . 




(4) 




(5) 
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The field tensors of the vector mesons and of the electromagnetic field are 
given by: 

W = d^uo v - d u uj» 

R/*" = d^p u - d u p^ (6) 

The variational principle gives rise to the Dirac equation: 

{-ia.V + V(r) + (3m*) ^ = e^, (7) 

where V(r) represents the vector potential: 

1 - r 3 

V(r) = g^uj (r) + g p T 3 p (r) + e A (r), (8) 

and S(t) is the scalar potential 

S(r)=g a a(r). (9) 

The effective mass is defined by the scalar potential as 

m * (r) =m + S(r). (10) 

The corresponding Klein-Gordon equations can be written as: 



(- A + mf)a = -g„W> 
(-A + m* 2 )^ = g^ v i) 
(-A + mf)p u = g p i>~i v Ti> 

-AA U = \ei>{l + r 3 ) 7 ^. 



(11) 



For the linear Lagrangian model, 

m* = m a ; m* = m„; m* = m p . (12) 



Thus, the meson masses do not exhibit a density dependence in the linear 
model. For an even-even nucleus with time-reversal symmetry, the spatial 
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components of the vector fields u), ps and A vanish. The Klein-Gordon equa- 
tions for the meson fields are then time-independent inhomogeneous equations 
with the nucleon densities as sources: 



(-A + m>(r) = -g a p s (r) 
(-A + m> (r) = g u p v {r) 
(-A + m 2 p )p (r) = g p p 3 (r) 
-AA (r) = ep c (r). 



where the source terms p s , p v , p 3 and p p to the Klein-Gordon equations (Eq. 
[T51) are the scalar, vector, isovector and charge densities, respectively, as de- 
fined by nucleon spinors: 

A _ 

p s = E^i A 

A 

pv= E A 

l= z N (14) 

P3 = E - E Ipn 

p=l n=l 
Z 

Pc= E ip P - 

p=l 



Here, the sums are taken over the valence nucleons only. The stationary state 
solutions ipi are obtained from the coupled system of Dirac (Eq. [7]) and Klein- 
Gordon equations (Eq. [T5|) self-consistently. 

The total ground-state energy of a spherical nucleus can be expressed as a 
functional of the baryon spinors {ipi} 

Erm f [ifi] = (®\H\$), (15) 
obtained using the Hamiltonian density 

/ HRMF{r) = r(r) + mp s {r) 

+^9aps(r)(r(r) + -g u p v (r)uj (r) 

+ ^9 P p3(r)p°{r) + ^ep p {r)A°(r) 
-i((V^(r)) 2 + m>°(r)) 2 
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-\{(yA\r)f (16) 

which is given in terms of the source densities (Eq. [T4"]) . The kinetic energy 
density is given by 

T {r)=Y,4{r){-i°X}A{r). (17) 
i=i 

Taking the variation of Eq. ffl5l) with respect to ipj, the stationary Dirac 
equation (Eq. [7]) with energy eigenvalues is obtained: 

h D ^i{r) = eii>i(r), (18) 
where 

h D = -iaV + (3m* + g w cu°(r) + g p r 3 p°(r) + e (1 ~ ^ A°(r). (19) 

Solving the Dirac equation (Eq. [18]) self-consistently, the ground-state $ of 
the nucleus is written as a Slater determinant of single-particle spinors ipi (i 
= 1,2,...,A). 

The linear Walecka model has been successful in attaining saturation of nuclear 
matter as a delicate balance between large fields due to a and u mesons. 
However, a proper description of the properties of finite nuclei was not possible 
until nonlinear self-couplings of the a meson were introduced [31]. 



4 The RMF Lagrangian models 

In this work, we have considered three successful RMF Lagrangian models 
in order to analyze the breathing-mode isoscalar GMR. In the following, we 
provide a brief description of the formulation of the various Lagrangian models 
considered. 

4-1 The nonlinear-a model 

The nonlinear sigma (NLcr) model is the standard Lagrangian model that is 
used most commonly for calculation of the ground-state properties of nuclei. 
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An important ingredient of the NLa model is the assumption of nonlinear 
scalar self-couplings of the form [M] : 



U NL = \g2<y' A + \g^\ (20) 



The parameters $2 and g% are the nonlinear couplings of the a-meson in the 
conventional a 3 + a 4 model. The effective Lagrangian for the NLa model 
becomes 

£e// = £o _ Unl- (21) 

The Klein-Gordon equations (Eq. [TT]) produce a density dependence of the 
meson masses as given by 

mf = ml + 5- 2 o- + g3& 2 

m* = m, (22) 
m* = m p . 



Only the cr-meson mass exhibits an implicit density dependence in the NLcr 
model. 

The scalar self-couplings have proved to be important for an appropriate de- 
scription of nuclear surface and have thus become indispensable. However, the 
negative quartic coupling in the model NLu has been a source of instability in 
nuclear matter at higher densities P^ll35f36j . It has been shown recently [37] 
that a scalar-vector interaction (SVI) model that comprises a combination of 
couplings of a and u mesons can dispense with the a 3 + a 4 terms in the RMF 
Lagrangian. The SVI model is able to provide an improved description of the 
ground-state binding energies and charge radii of nuclei. We shall discuss the 
basic features of the SVI model in this section below. 

The NLcr model is well-established and has shown to be a successful model 
for calculating ground-state properties of nuclei. The earliest parameter sets 
obtained were NL1 (K^ = 211 MeV) [21] and NL2 (K^ = 399 MeV) |38]. 
However, due to the large asymmetry energy J, these sets were not appropriate 
for nuclei away from the stability line. In remedying the problem of a large 
neutron skin [31], the force NL-SH [2H] was developed as one of the first 
successful parameter sets in the RMF theory, which also described nuclei away 
from the stability line. Due to a relatively larger value of of NL-SH, the 
force NL3 has been obtained [29] with a view to getting a description of the 
breathing-mode GMR in a physically acceptable region. 
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Table 1 

The nuclear matter properties of various parameter sets in the NL<r model. 



Sets 


Koo (MeV) 


E/A (MeV) 


m* 


p (mi" 3 ) 


J (MeV) 


NL1 


211 


-16.42 


0.573 


0.152 


43.5 


NL3 


272 


-16.25 


0.595 


0.148 


37.4 


NL-SH 


355 


-16.32 


0.597 


0.146 


36.1 


NL2 


399 


-17.02 


0.667 


0.146 


45.1 



In this work, we consider the forces NL1, NL3, NL-SH and NL2, whose order 
reflects an increasing value of K^. Nuclear matter properties of these param- 
eter sets are provided in Table 1. It may be noted that even though the forces 
NL1 and NL2 are quite different in some of the properties such as incom- 
pressibility and asymmetry energy, these are able to reproduce ground-state 
properties of nuclei along the stability line quite well. Here, we have included 
NL2 in order to extend the systematic behaviour of the NLcx model into the 
region of high incompressibility. 



4-2 The a-u coupling model - SIGO 



The a-uj coupling model (SIGO) has recently been introduced by Haidari and 
Sharma [40j with a view to bring about an improvement in the ground-state 
properties of nuclei. A coupling between a and u mesons of the form 



U„ 



1 1 
:9auO- to 



(23) 



was introduced in addition to the usual NLa scalar potential of the form a 3 
+ a 4 . The effective Lagrangian for the model SIGO becomes 

C eff = C -U NL + U acJ . (24) 



The meson mass terms in the corresponding Klein-Gordon equations (Eq. [IT!) 
are then given by: 

mf = ml + g 2 a + g 3 a 2 - g auj u 2 

m *J = m l + 9<ru(r 2 (25) 
m* = m p . 
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Thus, in the SIGO model, both the a and 00 meson masses exhibit an implicit 
density dependence. 

It has been shown [ID] that the parameter set SIG-OM obtained in the model 
SIGO is able to provide an excellent description of the ground-state properties 
such as binding energies and charge radii of nuclei along the stability line as 
well as far away from it. Especially, a significant improvement in the binding 
energy of nuclei at the magic numbers has been achieved. Consequently, an 
excellent description of the total binding energy of Sn isotopes all over the 
range of the shell from 100 Sn (N = 50) to 132 Sn (N = 82) was obtained. 
Charge radii of nuclei, especially of Pb isotopes, which are over estimated 
significantly by NL3, were well reproduced. The ensuing EOS of the nuclear 
matter with SIG-OM has been shown to be much softer than that given by 
the parameter sets of the NLcx model [40J. 

In order to investigate the behaviour of the model SIGO for the breathing- 
mode GMR, we have constructed a few parameter sets by extending the range 
of below and above that of the set SIG-OM (K^ = 265 MeV). The 
parameter sets SIGO-a, SIGO-b, SIGO-c and SIGO-d have been obtained 
with ranging between 241 - 283 MeV. The nuclear matter properties of 
the parameter sets of the model SIGO are given in Table 2. It may be noted 
that the sets SIGO-a, SIGO-b, SIGO-c and SIGO-d corresponding to different 
values of are able to reproduce the binding energy and charge radii of key 
nuclei from 16 to 208 Pb well. Qualitatively, the set SIG-OM is considered as 
the best amongst all the sets of the SIGO model provided in Table 2. 



Table 2 

The nuclear matter properties of the a-co coupling (SIGO) sets. 



Sets 


Koo (MeV) 


E/A (MeV) 


m 


Po (fm" 3 ) 


J (MeV) 


SIGO-a 


241.1 


-15.90 


0.621 


0.149 


33.0 


SIGO-b 


248.4 


-16.00 


0.622 


0.150 


33.6 


SIG-OM 


265.2 


-16.30 


0.622 


0.149 


37.0 


SIGO-c 


272.8 


-16.33 


0.623 


0.149 


37.7 


SIGO-d 


282.5 


-16.34 


0.620 


0.148 


37.5 
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4-3 The scalar-vector interaction model - SVI 



The scalar-vector interaction (SVI) model has recently been developed by 
the author |37j with a view to remove nonlinearities of the mesonic fields 
in the RMF theory. Self-interactions of a field have been cited as a source 
of instability in the nuclear matter at higher densities [36]. It has now been 
demonstrated [37] that by a suitable combination of couplings of a and uj fields, 
the self-interactions of the mesonic fields can be dispensed with. It is expected 
that the introduction of SVI would be consistent with a linear realization of 
the chiral symmetry in a suitable scheme [H]. 

The SVI model consists of the meson-meson interactions between a and uj 
mesons of the form 

Umra = \g^uj 2 + ]-g 5 cr 2 uj 2 , (26) 



where g± and g$ are the respective coupling constants for interactions be- 
tween a and uj mesons. The effective Lagrangian without self-interactions of 
the bosonic fields thus becomes 

£eff = ^0 + Umm- (27) 



The corresponding Klein-Gordon equations (Eq. [TTj) have the effective meson 
masses as given by: 

m* 2 = ml- g 4 oj 2 /(2a) - g 5 uj 2 

m *J = m l + 9aO- + g 5 a 2 (28) 
m* = m p . 



These equations exhibit an implicit density dependence of the a and u meson 
masses. Both the a and uj meson masses provide a density dependence that 
is different from that for the SIGO model (Eq. [251) . Thus, the density depen- 
dence of a and uj meson masses, as represented by Eqs. (1221 . (|25|) and (|28|) for 
the three Lagrangian models NLa, SIGO and SVI, respectively, are different 
from one another. A comprehensive discussion of the differences in the den- 
sity dependences in various Lagrangian models and their influence on nuclear 
properties shall be presented elsewhere [13]. It shall be instructive to see as 
to how these density dependences in the high-density regime would compare 
with the predictions of the Brown- Rho scaling [44] . 
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Table 3 

The nuclear matter properties of the scalar-vector interaction (SVI) sets. 



Sets 


K (MeVl 


El A fMeVl 


* 

ID 
1 1 1 


f)r\ ( fm ^\ 
fJQ yim J 


j (MpV) 


SVI-a 


243.5 


-16.30 


0.615 


0.150 


37.4 


SVI-b 


253.1 


-16.24 


0.612 


0.150 


36.6 


SVI-1 


263.9 


-16.30 


0.616 


0.149 


37.6 


SVI-2 


271.5 


-16.31 


0.621 


0.149 


37.0 


SVI-c 


288.1 


-16.25 


0.617 


0.148 


36.4 



The parameters sets SVI-1 (K^ = 264 MeV) and SVI-2 (K^ = 272 MeV) 
were obtained in Ref. [37]. These sets have been shown to provide an excellent 
description of the ground-state energies of nuclei over a large range of the 
periodic table. This includes nuclei along the stability line as well as far away 
from it. Description of binding energies and charge radii of nuclei that is 
achieved with the SVI sets is significantly better than that with the NLcx set 
NL3. Having established the basis of the SVI model, we have constructed a few 
more parameter sets having in the range ~ 240 — 290 MeV for a systematic 
analysis of the breathing-mode GMR. The nuclear matter properties of the 
parameter sets SVI-a, SVI-b, and SVI-c thus constructed in addition to those 
of SVI-1 and SVI-2 are listed in Table 3. It is worth mentioning that the 
auxiliary sets SVI-a, SVI-b, SVI-c are also able to provide a reasonably good 
description of the ground-state properties of nuclei. 



5 The generator coordinate method and the breathing-mode GMR 

The generator coordinate method (GCM) is a powerful tool to study ground 
and excited states in atoms and nuclei. The GCM has been applied to study ef- 
fects of correlations on the ground-state properties of nuclei using the density- 
dependent Skyrme approach [35]. In the present work, we employ the GCM in 
the RMF theory to explore the excited state of the isoscalar GMR in atomic 
nuclei. 

The GCM is based upon a trial A-particle wavefunction ^gcm in the form of 
a linear combination of 

*gcm(*i, r 2, r A ) = / J r (q) ( ^(r 1 , r 2 , ...,r A ;q)dq, (29) 
where the generator function = $(r 1; r 2 , r A ; q) is a Slater determinant 
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constituted from single-particle spinors if)i(r, q),(i = 1, 2, A) as a func- 
tion of the generator coordinate q. The "weight" or the "generator function" 
J-{q) is determined from variation of the total energy of the system 



E[F\ = ^ GCM| , '" GCW (30) 



(^GCM 






^gcm) 



with respect to T(q). This leads to the Hill- Wheeler equation for the weight 
function: 

/ Mr; q, q') - E^(r; q, q')]F(q')dq' = 0, (31) 



where 

n{r ] q 1 q') = ^{q)\H\^{q')) (32) 

is the overlap energy- density kernel of the Hamiltonian H associated with the 
RMF Lagrangian, and 

A/Xr;g,<zO = <*(?) W)> (33) 

is the overlap norm kernel. The overlap energy-density kernel takes the form 
given in Eq. (IIEl) : 

H RM F(r; q, q')=r(r; q, q) + Mp s (r; q, q) 

+ ^9aPs{r\ q, q')a(r; q, q') + -g u p v (r; q, q')uJ°(r; q, q') 

+^9 P P3(r; q, q')p°(r; q, q') + ^ep p (r; q, q')A°(r) 
+-((Vo"(r; q, q')f + m^a 2 (r; q, q') 
- l -{{Vu\r-q,q')f + ml{ui\r-q,q>)Y 
-l((Vp (r;q,q')) 2 + m%p°(r;q,q>)) 2 
-i((VA°(r;g,g')) 2 

+ other terms. (34) 

Here, the 'other terms' refer to the meson-meson interaction terms of various 
Lagrangian models considered in Section HI The kinetic energy density can be 
written in terms of the spinors {^(r; q)} as 

r(r;q,q') = X^Vfe g){-iaV}V;(r; g') (35) 
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The source densities appearing in Eq. fl34|) can be written as 



p s (r;q,q') = ^N-^^Y-q^^r-q'), 
Pv(r;q,q') = ^N^tpjfaqtyjfaq'), 
P3(r;q,q') = J2 N ii Vi( r ; <?)^j(r; q'), 

Pp (r;q,q>) = ^N^j(Y;q)^^-^(r,q'). (36) 
The overlap matrix elements (q, q') are calculated using the spinors as 
Na(q, q') = J rf 3 r ^(r; g)%(r; q') (37) 



The determinant of Nij(q, q') provides the overlap kernel (Eq. 1331) : 

Ar(q,q')=det{N tJ (q,q > )}. (38) 



On obtaining the energy and norm overlap kernels using Eqs. (I3"2j) and 
the Hill- Wheeler equation (13"T|) is solved. The solution provides the nuclear 
ground and the excited states. 

In order to obtain the isoscalar GMR energy of nuclei, we perform constrained 
RMF calculations by solving the Dirac equation 

(h D - qQ)ipi(x) = eiipi(x) (39) 



with the constraint operator Q = r 2 . For the case of the isoscalar GMR, the 
Lagrange multiplier q is associated with the value of the nuclear root-mean- 
square (rms) radius 

R=(Q) = {jfr 2 p v (r-q)d 3 r} 1/2 (40) 



where p v (r; q) is the baryon density determined by the solution {^i(r; q)} for 
a given value of the generator coordinate q. 

Constrained RMF calculations have been performed in coordinate space for 
a range of the Lagrange multiplier q. Constructing the generator Slater de- 
terminant $(g), the integral kernels, Eq. (l32j) and Eq. (|33|) are obtained. On 
solving the Hill- Wheeler equation, the first excited state gives the energy of 
the isoscalar GMR. Details of the GCM method applied to the RMF theory 
have been provided in Ref. 
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6 Results and discussion 



As noted earlier in Section HJ the Lagrangian sets SIG-OM in the SIGO model 
[iO] and SVI-1 and SVI-2 in the SVI model [37] developed recently have been 
shown to bring about a significant improvement in the ground-state properties 
of nuclei as compared to the parameter set NL3 of the NLcr model. With the 
advent of the alternative Lagrangian models SIGO and SVI, we have now 
explored the response of various Lagrangian models to the breathing-mode 
isoscalar GMR. 

As alluded to in Ref. [ID], the set SIG-OM yields GMR energies which are 
larger than those of NL3, though the incompressibility iV» of SIG-OM is 
smaller than that of NL3. It seems paradoxical, for according to conventional 
wisdom, a force with a higher value of should yield a higher value of the 
GMR energy and vice versa. This strange feature has prompted us to investi- 
gate the GMR response within the framework of the RMF theory employing 
various Lagrangian models at hand. The isoscalar GMR energies have been 
calculated for a few key nuclei using the GCM approach. 




261 266 271 

(MeV) 

Fig. 1. The breathing-mode GMR energies for nuclei 208 Pb, 120 Sn and 90 Zr obtained 
with the generator coordinate method using the NLcr Lagrangian set NL3, the SVI 
Lagrangian sets SVI-1 and SVI-2 and the SIGO Lagrangian set SIG-OM. The sets 
NL3 and SVI-2 have nearly the same value of K^. 
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6.1 The GMR energies: a paradoxical picture 

In order to visualize the response of various Lagrangian models, we show the 
GMR energies obtained by employing the GCM approach in Fig. [T] The GMR 
energies for three key nuclei of 208 Pb, 120 Sn and 90 Zr have been obtained using 
the parameter set NL3 of the NLcx model, with the sets SVI-1 and SVI-2 of the 
SVI model [37] and with the set SIG-OM of the SIGO model [ffi]. Contrary to 
the usual increase of the GMR energy with Koo, an unusual pattern displaying 
rather paradoxical behaviour is seen. The GMR energies obtained with SIG- 
OM are bigger than those of NL3 as mentioned above. These are also bigger 
than those of SVI-2. Moreover, SVI-2 values are bigger than those of NL3, 
though SVI-2 and NL3 have nearly same value of K^. This behaviour implies 
that some finite-size contribution(s) to the incompressibility must be different 
in these Lagrangian models. Such a behaviour was noted in relativistic RPA 
calculations in Ref. [IB]. It was shown that the force NL-SH of NLa model 
and TM1 with quartic to coupling [47] . both having very different values of 
Koo, were shown to provide comparable values of the GMR energy for 208 Pb. 
In order to throw light upon the seemingly paradoxical behaviour of the GMR 
energies and with a view to gauge the response of the finite-size effects to 
the nuclear incompressibility of finite nuclei, we have first analyzed the GMR 
energies within each Lagrangian model separately. 



6.2 GMR energies with the NLo model 

In the NLcr model, we have employed the standard parameter sets NL1, NL3, 
NL-SH and NL2, which are in an increasing sequence of K^. These span a large 
range of ~ 210-400 MeV. The isoscalar (T=0) GMR energies calculated 
for nuclei 208 Pb, 120 Sn and 90 Zr using these parameter sets are shown in Fig. 
[2j The GMR energies are shown directly as a function of in Fig. EJja). Due 
to obvious dependence of E GMR on Ka, the curves exhibit a quadratic like 
dependence on K^. Only for 90 Zr does the energy seems to fall slightly out 
of the trend especially for NL2. This may be due to an extremely large value 
of Koo ~ 400 MeV for NL2, whereby some anharmonicity may creep in for 
90 Zr, which is relatively light as compared to 208 Pb. It is interesting to note 
that the curve for each nucleus and especially that for the heavy nucleus 208 Pb 
shows a monotonic dependence on K^, though each of the parameter sets was 
constructed under different circumstances. Some nuclear matter properties of 
these sets such as the asymmetry energy are known to be rather different 
especially that of NL1 and NL2 with a value of J ~44 MeV. 

The GMR energies are shown as a function of \[K~ a Z ) in Fig. EJ^b) in order to 
eliminate the apparent quadratic dependence. The curve for 208 Pb shows a 
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Fig. 2. The breathing- mode isoscalar GMR energies for 08 Pb, °Sn and 90 Zr ob- 
tained with the standard parameter sets NL1, NL3, NL-SH and NL2 of the NLcr 
model in an increasing sequence of K^. The energies Eqmr are displayed as a 
function of (a) ifoo and (b) \fK^. 



slight improvement towards a linear behaviour. For Sn and Zr a similar 
feature may not be visible due to interplay of stronger finite-size effects such 
as the surface which becomes important for medium heavy and lighter nuclei. 



6.3 The GMR energies with the SVI model 



We have calculated the breathing-mode GMR energies within the scalar- vector 
Lagrangian model SVI using the GCM with the parameter sets SVI-a, SVI-b, 
SVI-1, SVI-2 and SVI-c, which are in an increasing order of K^. These sets 
encompass a region of ~ 240 — 290 MeV. The results are shown in Fig. [3J 
The curves show a monotonous behaviour with K^. All the parameter sets of 
SVI fall in an orderly pattern and show a dependence on similar to that 
seen in Fig. [2j However, as compared to the values in Fig. [2J the GMR energies 
with SVI are higher than the corresponding values with the NLcr model. 
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Fig. 3. The breathing-mode isoscalar GMR energies for 08 Pb, °Sn and 90 Zr 
obtained with the parameter sets SVI-a, SVI-b, SVI-1, SVI-2, and SVI-c of the 
scalar- vector interaction model S VI of Ref. [37] . 

6.4 The GMR energies with the SIGO model 



We have then investigated the response of the a-u model SIGO to the breathing- 
mode GMR. GCM calculations have been performed with the parameter sets 
SIGO-a, SIGO-b, SIG-OM, SIGO-c and SIGO-d. The results are shown in Fig. 
[H All the SIGO sets display an orderly pattern such as that seen in Figs. [2] 
and [3] with the NLcr and SVI models, respectively. Comparatively, the SIGO 
values are higher than the corresponding SVI values. 



6.5 A comparative picture 



The results obtained with the three Lagrangian models are compared in Fig. 
Focusing on the region of physically acceptable values of K^, the figure shows 
the GMR energies for the three Lagrangian models within the range of ~ 
230—320 MeV. For a given value of K^,, the model NLcr gives the lowest values, 
whereas the model SIGO delivers the largest values for the breathing-mode 
GMR energy. Thus, the models NLcr, SVI and SIGO, respectively, produce 
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Fig. 4. The GMR energies for 208 Pb, 120 Sn and 90 Zr obtained with the parameter 
sets SIGO-a, SIGO-b, SIG-OM, SIGO-c and SIGO-d of the a-u model SIGO of Ref. 

m 

GMR energies in an increasing order for any given nucleus. The results in 
this figure demonstrate that the GMR energies depend strongly upon the 
Lagrangian model employed and that each model predicts a different value 
of the GMR energy for a given K^. This implies that the GMR energy is 
not a simple function of only the incompressibility of the infinite nuclear 
matter. Different values of the GMR energy for a given value of indicate 
different contribution of finite-size effect (s) to the incompressibility Ka of a 
nucleus. We discuss below the relative magnitudes of various contributions 
to Ka and attempt to discern the factors responsible for multiplicity of the 
GMR energies for a given value of and thus the reason for producing the 
paradoxical behaviour seen in Fig. [IJ 

The multiplicity of the GMR energies for a given in the RMF theory is in 
contrast to a single value that is obtained usually within the density-dependent 
Skyrme theory [14] . This can be attributed to a difference in the density de- 
pendence of the nuclear interaction amongst the various Lagrangian models. 
In comparison, the density dependence within the non-relativistic Skyrme ap- 
proach is provided by the standard Skyrme density functional used universally 
- which is fixed at the outset. Should an alternative form of the density func- 
tional in the Skyrme approach arise, a situation akin to the RMF theory would 
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emerge. Microscopic differences amongst various Lagrangian models, which 
lead to a phenomenal difference in the GMR energies, are being investigated 
separately [33] . 




i i i i i i i i 

230 270 310 



(MeV) 

Fig. 5. The GMR energies for 208 Pb, 120 Sn and 90 Zr for the three Lagrangian models 
in the physically acceptable region of compared with the experimental data. The 
arrows show intersection of the experimental data (dashed lines) with the theoretical 
curves obtained with the NLcr model. The results with sets having nearly the same 
value of i^oo ~ 272 MeV in the three Lagrangian models are enclosed by rectangular 
boxes. 

6.6 Theoretical extraction of - a conundrum 

Given the emerging situation in the RMF theory, an intersection of the ex- 
perimental values (shown by the dashed horizontal lines in Fig. \5§ with the 
theoretical curves would yield a different value of depending upon the 
Lagrangian model used. The experimental values shown by the dashed lines 
cross the curves at different locations in K^. Staying with the pivotal case of 
208 Pb in Fig. [5](a), each Lagrangian model yields value which is different 
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from the others. inferred from the datum of 208 Pb with the model NLcr 
is ~ 310 MeV, whereas a value of ~ 290 MeV and ~ 260 MeV would be 
determined with the models SVI and SIGO, respectively. 

Using the datum on 120 Sn in Fig.[5](b), the value of inferred is ~ 300 MeV, 
~ 275 MeV and ~ 245 MeV with the NLcr, SVI and SIGO models, respectively. 
Similarly, using the datum on 90 Zr in Fig. 0(c), one would conclude a value 
of as ~ 280 MeV, ~ 265 MeV and ~ 240 MeV, respectively. The shift 
in the value of inferred from a given model using the experimental data 
on 208 Pb, 120 Sn and 90 Zr is illustrated by the arrows for the case of the NLa 
model. A similar shift is evident in the figure (not shown by arrows) for the 
models SVI and SIGO as well. 

Staying with the interpolation method, it is equally illustrative to conclude 
that any of the Lagrangian models is not able to match the experimental 
data on all the three nuclei simultaneously with a given value of K^. These 
observations bring into question the hitherto used theoretical approach based 
upon interpolation amongst forces. Thus, the approach of determining a la 
Blaizot [1] becomes unsustainable in the RMF theory. It is therefore essential 
to 'calibrate' the density dependence of the nuclear interaction in the RMF 
theory before such an approach could be applicable. Arguably, the form of 
the density dependence opens up another dimension in the landscape of the 
breathing-mode GMR and the incompressibility of nuclear matter. 

The shift of the arrows from the right towards the left in going from 208 Pb 
to 90 Zr (i.e., a decrease in the value of so inferred) needs analysis. It is 
understood that the interior (or bulk) of the nucleus 208 Pb has a density close 
to that of the saturation density of nuclear matter and hence is synonymous 
with the nuclear matter for practical purposes. It has a nuclear surface that 
is smaller as compared to that in 120 Sn and 90 Zr. Thus, in going from 208 Pb 
to 90 Zr, the region of the nuclear surface increases. The effect of compression 
of a correspondingly large nuclear surface needs to be accounted for. The 
continuous shift in the arrows (and hence value inferred) from 208 Pb to 
90 Zr points towards the increased role played by an increasingly larger surface 
region. 



6. 7 Disentangling the finite-size effects and the surface incompressibility 

A shift in the value of extracted from the experimental GMR energies of 
208 Pb, 120 Sn and 90 Zr, as shown in Fig. [5], renders the approach of interpolation 
amongst forces as futile. It is, therefore, pertinent to revisit the expansion of 
Ka to dissect various components. As noted earlier, the use of Eq. (j2J) does not 
allow an unambiguous determination of primarily due to the correlation of 
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the Coulomb term to superimposed on to a much smaller magnitude of the 
former. Here, we adopt a slightly different approach. Instead of allowing the 
aforesaid correlation to jeopardise the extraction of K^, we let the coefficient 
Kcoui fit as an independent variable. 

For a comparative analysis of the GMR energies, we have selected the pa- 
rameters sets having the same value of in the three different Lagrangian 
models. The sets with an incompressibility close to ~ 272 MeV are en- 
closed by the rectangular boxes in Fig. [51 This corresponds to the sets NL3 
(K^ = 271.8 MeV), SVI-2 = 271 A MeV) and SIGO-c (K^ = 272.8 

MeV) in the Lagrangian models NLo~, SVI and SIGO, respectively. 

Table 4 

The GMR energies for 208 Pb, 90 Zr, 48 Ca, and 40 Ca obtained with the parameter 
sets NL3, SVI-2 and SIGO-c having the same value of ~ 272 MeV in the three 
Lagrangian models. 



Nucleus 


NL3 


SVI-2 


SIGO-c 


208p b 


13.0 


13.5 


14.2 


90 & 


16.9 


17.5 


18.4 


48 Ca 


18.9 


19.7 


20.5 


40 Ca 


19.6 


20.3 


21.1 



We have selected the key nuclei of 208 Pb, 90 Zr, 48 Ca, and 40 Ca spanning a 
broad range of masses for our analysis. 48 Ca and 40 Ca are included with a 
view to augment the range of variation with A and to reinforce the asymmetry 
component, and thus to facilitate a reasonable dissection of the components in 
the analysis. As observed in Refs. [T2ll21f48] , Sn isotopes present an anomalous 
behaviour and are at present problematic theoretically. Sn nuclei being open- 
shell exhibit a significant BCS neutron superfluidity in the ground state thus 
affecting the Fermi surface significantly. These may require further theoretical 
considerations than employed so far. With the problem of Sn nuclei not yet 
understood, we have not included 120 Sn isotope in the dissection of finite-size 
contributions. 

The GMR energies for 208 Pb, 90 Zr, 48 Ca, and 40 Ca obtained with the GCM us- 
ing the parameter sets NL3, SVI-2 and SIGO-c of the three Lagrangian models 
are given in Table 4. As seen earlier, NL3, SVI-2 and SIGO-c, respectively, 
produce GMR energies in an increasing order for all the nuclei considered. The 
energies with NL3 are the lowest, whereas SIGO-c provides the largest values 
amongst the three Lagrangian models. A preliminary consideration suggests 
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that Koo being the same for these parameter sets, the finite-size contributions 
to Eq. (j2J) must be largest with NL3 and smallest with SIGO-c in order to 
have the noted difference in the GMR energies. 

In order to discern which finite-size effect (s) play a crucial role in influencing 
the GMR energies such as given in Table 4 and portrayed in Figs. we 
have resorted to the use of Eq. The objective is to see whether a fit to 
Eq. (|2D is possible by relinquishing the aforesaid correlation of Kc ou i to K^. 
Using the GMR energies given in Table 4 for the parameters sets of the three 
Lagrangian models, we have performed a ^-minimization procedure to fit Eq. 
d2J) by treating K sur f, K asym , K Cou i and K curv as independent variables. The 
corresponding Ka is evaluated using Eq. (TTJ. 

Table 5 

The coefficients of expansion (in MeV) of Ka for the parameter sets of the three 
Lagrangian models with the same value of ~ 272 MeV. The last column shows 
the ratio K sur j / Kqo obtained. 



Model 




^surf 


K 


Kcoul 


K 

1 1 curv 


Ksurf 1 


NL3 


271.8 


-537 ± 15 


-389 ± 12 


-6.9 ±0.2 


117 ±43 


-1.98 


SVI-2 


271.4 


-455 ± 15 


-345 ± 12 


-6.8 ±0.2 


-52 ± 45 


-1.67 


SIGO-c 


272.8 


-273 ± 17 


-395 ± 13 


-7.1 ±0.3 


-622 ± 50 


-1.00 



For the purpose of this work, i.e. to investigate the importance of finite-size 
effects, we fix the value of to the theoretical one and allow a fit of the GMR 
energies only to four free parameters, viz., K sur f, K asym , Kc ou i and K curv . It 
is noted that by doing so it is possible to fit Eq. (j2J) to the GMR energies well. 
The results of the fits of the GMR energies for the three Lagrangian models 
are given in Table 5. The resulting coefficient Kcoui turns out be similar in 
value from NL3 to SIGO-c. It is not expected that the Coulomb term would 
be very different from one Lagrangian model to the other. These values are 
very close to those which have been obtained using other approaches such as 
the density-dependent Skyrme forces [20J. 

It can be seen from Table 5 that by fixing to the desired value, the 
surface term can be determined with a considerable precision. The striking 
feature that emerges from the minimizations is the significant difference in 
the value of the surface incompressibility. Here K sur f is obtained as —537 
MeV for NL3, -455 MeV for SVI-2 and -273 MeV for SIGO-c. The set NL3 
provides the largest value of K sur f, whereas SIGO-c yields the lowest value. 
Inevitably, a large difference in the value of K sur f for the three Lagrangian 
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models is responsible for the differences in the GMR energies as observed in 
Figs. [5] and [61 This is also reflected by the significant differences in the ratio 
K sur f I 'Koq shown in the last column. It varies from —1.98 for NL3, —1.67 for 
SVI-2 to -1.00 for SIGO-c. The ratio with NL3 and SVI-2 is substantially 
different from the ratio of ~ — 1 that is obtained from calculations of semi- 
infinite nuclear matter using the scaling model [42J. On the other hand, a ratio 
of ~ —1 is attained with SIGO-c. 

In comparison, the asymmetry coefficient K asym shows only a limited variation 
from ~ —345 MeV to ~ —395 MeV (Table 5). These values are not far from the 
value of —550 ± 100 MeV obtained in a recent fit of the experimental GMR 
energies of Sn isotopes [12]. A value of —500 ± 50 MeV has been obtained 
in a recent analysis of the breathing-mode GMR with the density-dependent 
Skyrme forces [20J. 
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Fig. 6. Contributions of various terms in the expansion Ka = T,ciKi using the 
coefficients K{ from Table 5 for the parameter set SVI-2. 

For the sake of illustration of the relative magnitudes of various contributions 
to Ka in Eq. (|2l) , we show the size of various terms using the coefficients 
obtained for the parameter set SVI-2 in Fig. [61 Here we have chosen the 
set SVI-2 for the sake of illustration. The picture with the other parameters 
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sets is similar to it with major differences being in the surface contribution. 
The curvature term would show differences as maifest in Table 5 especially 
for SIGO-c. The variation of the Coulomb and the asymmetry terms with 
the mass number A is relatively small. Additionally, the Coulomb and the 
asymmetry contributions are an order of magnitude smaller than the surface 
term. 

The surface term provides the largest contribution to Ka amongst the finite- 
size effects. It also shows a significant variation with A. Consequently, the 
surface incompressibility has the ability to modulate the GMR energies of 
medium mass and lighter nuclei in a significant measure. 
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Fig. 7. The surface contribution to Ka for the parameter sets NL3, SVI-2 and 
SIGO-c in the three Lagrangian models having the same value of ~ 272 MeV. 

The relative magnitudes of the surface term for the parameter sets of the three 
Lagrangian models having the same value of ~ 272 MeV are compared 
in Fig. [71 The surface term exhibits a large variation with mass A. Significant 
differences in the surface term amongst the various Lagrangian models are seen 
readily in the figure. Inevitably, these differences give rise to the differences in 
the breathing-mode energies obtained with various Lagrangian models. 

Given the relatively smaller magnitude of the Coulomb and the asymmetry 
terms as seen in Fig. [6j magnitudes of these terms would not show any signif- 
icant difference from one Lagrangian model to the others. In conclusion, it is 
the surface term that plays a decisive role in influencing the GMR energies of 
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nuclei amongst various Lagrangian models significantly. Whilst the curvature 
term is not very small especially with SIG-c, it may also be influencing the 
GMR energies to an extent. 



6.8 Density dependence of the meson masses 

The three Lagrangian models employed in this work are successful in reproduc- 
ing the ground-state properties of nuclei. Whilst most of the nuclear matter 
properties e.g., saturation density, saturation binding energy, effective mass 
and asymmetry energy are not very different from one model to the other, the 
density dependence of the a and uj meson masses does show different forms in 
various models as seen in Section HI This difference in the density dependence 
especially in the sub- saturation densities would influence the surface proper- 
ties. Consequently, the form of the density dependence of meson masses affects 
the GMR response implicitly. 

The density dependence of meson masses and by implication that of the nu- 
clear interaction in the RMF theory remains an open and challenging problem. 
The introduction of the nonlinear a 3 + a 4 terms in the RMF Lagrangian has 
led to successful results on ground-state properties of nuclei. It has, however, 
a major drawback in that the resulting EOS of the NLcx model is ubiquitously 
stiff. The standard NLcx sets produce a maximum mass of neutron stars in 
the vicinity of 2.7-3.0 solar masses. This is well above the spectrum of the 
observed neutron star masses. Thus, the NLcr model is not deemed as suitable 
for neutron star structure. In view of this, development of other Lagrangian 
models which could describe finite nuclei as well as nuclear matter at higher 
densities becomes desirable. 

The recent development of the SVI [37j and SIGO [40J models provide a much 
improved description of nuclear properties vis-a-vis the NLa model. Notably, 
the SIGO model provides an excellent description of nuclei at the shell closures 
and also for nuclei in the extreme regions of the periodic table [40] . It also 
delivers an EOS of nuclear matter which is softer than the NLa counterparts. 
On the other hand, the Lagrangian model SVI is able to describe the ground- 
state properties of nuclei without self-interactions of the mesonic fields. It has 
been shown that introduction of terms of scalar-vector type such as those 
in SVI would help alleviate the onset of instabilities in nuclear matter |36j . 
Arguably, the SVI model would also be conducive to a linear realization of 
chiral symmetry. 

Notwithstanding the above, the RMF theory offers alternative possibilities of 
its density functional to be applicable to nuclei and nuclear matter. A natural 
question that arises is: what is the density dependence of the meson masses 
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that is preferred by nature? At present we do not have an answer to this 
question. It is difficult to pinpoint as to which Lagrangian model is preferred 
by experimental data. A comprehensive investigation of the RMF theory with 
all its attendant problems, virtues and theoretical constraints juxtaposed to a 
host of experimental data would be necessary in order to be able to unravel 
the form of the density dependence required by nature. 

6. 9 The established parameter sets and confrontation with GMR energies 

In view of the desirability of discerning the density dependence in the RMF 
theory, it is instructive to examine the results of the GCM calculations on 
breathing-mode GMR using the standard parameter sets vis-a-vis the available 
experimental data. The parameter sets which come into such a consideration 
are the ones which are in the region of physically acceptable values of K^. 
This includes the sets NL1 and NL3 of the NLcx model, sets SVI-1 and SVI-2 of 
the SVI model and SIG-OM of the SIGO model. With the exception of NL1, 
all the other sets have been shown to describe binding energies and charge 
radii of nuclei along as well as far away from the stability line. 

The GMR energies of the key nuclei 208 Pb, 120 Sn and 90 Zr obtained with 
the GCM approach using the above parameter sets are shown in Table 6. 
Experimental data on these nuclei are well established and are shown for a 
comparison. The parameter sets are listed in an ascending order of the GMR 
energies. 

Table 6 

The breathing mode GMR energies obtained with constrained GCM calculations 
using the established parameter sets. The experimental data [10|15j are shown for 
comparison. The value (in MeV) for the sets is shown in parentheses. 



Sets 


NL1 


NL3 


SVI-1 


SVI-2 


SIG-OM 


exp. 




(211) 


(272) 


(264) 


(272) 


(265) 




208pk 


11.0 


13.0 


13.3 


13.5 


14.1 


13.96 ± 0.28 


120 Sn 


12.7 


15.0 


15.2 


15.4 


16.2 


15.52 ±0.15 


90 Zr 


14.1 


16.9 


17.2 


17.5 


18.2 


17.81 ±0.30 



Without resorting to a disentanglement of Ka into its respective components, 
we compare the GMR energies due to various parameter sets directly with 
the experimental values. The GMR energies with NL1 are ~ 3-4 MeV smaller 
than the experimental values. NL3, on the other hand, provides GMR energies 
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which are ~ 1 MeV smaller than the experimental ones for 208 Pb and 90 Zr. 
For 120 Sn, the difference of the NL3 GMR energy with the experimental value 
is ~ 0.5 MeV. 

In comparison, the SVI-1 GMR energies are slightly bigger than those of NL3 
and hence are closer to the experimental values than NL3. However, even 
with SVI-1, the difference with the experimental data, especially for 208 Pb 
and 90 Zr is not insignificant. On the other hand, the SVI-2 GMR energies for 
all the three nuclei come much closer to the experimental bounds. As is the 
case with SVI-1, the GMR energies for 208 Pb and 90 Zr with SVI-2 are on the 
lower side of the experimental values, whereas for 120 Sn, the SVI-2 value is 
in good agreement with the experimental datum. Overall, the SVI-2 results 
show a reasonable good agreement with the experimental data for the three 
nuclei. The SIG-OM value, in comparison, comes closer to the experimental 
one for 208 Pb. However, SIG-OM overestimates the experimental values for 
120 Sn by ~ 0.7 MeV and slightly for 90 Zr. Thus, the behaviour of Sn isotope 
remains anomalous in the GCM approach, as observed in nonrelativistic and 
relativistic RPA approach [20~f2T] . 




From the point of view of the ground-state binding energies and charge radii 
of nuclei and a good description of data along the stability line and away 
from it, in conjunction with the breathing- mode GMR energies, the parameter 
set SVI-2 of the SVI model without self-interactions comes closest to the 
empirical data. It has to be seen in further exploration of other aspects of 
nuclear structure and properties as to how close to the data the set SVI-2 will 
appear. It should be pointed out that this conclusion drawn here is subject 
to the GCM treatment of the breathing mode GMR. The constrained GCM 
approach delivers GMR energies slightly smaller than those obtained within 
the RPA approach. It is worth investigation as to how different the properties 
of the breathing-mode GMR would appear within the framework of relativistic 
RPA using various Lagrangian models. 



7 Summary and conclusions 

The breathing-mode GMR is investigated within the framework of the RMF 
theory using the generator coordinate method. The GMR energy has been 
calculated for a few key nuclei of 208 Pb, 120 Sn, 90 Zr, 48 Ca and 40 Ca covering a 
broad range of atomic mass. Using the Lagrangian set NL3 of the NLcx model, 
sets SVI-1 and SVI-2 of the SVI model and set SIG-OM of the SIGO model, a 
paradoxical behaviour of the GMR energies has been found. Contrary to the 
received wisdom, parameter sets with a higher value of are observed to 
yield a lower value of the GMR energy as compared to a higher GMR energy 
attained with sets having a lower value of K^. 
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In order to resolve the apparent paradox of the GMR energies in the RMF 
theory, the GMR response has been investigated using the three different La- 
grangian models viz., the NLcr model, the SVI model and the SIGO model. 
Employing the constrained GCM approach, GMR energies have been calcu- 
lated within each Lagrangian model using parameter sets which encompass a 
broad range of K^. It is shown that each Lagrangian model exhibits a GMR re- 
sponse that is different from the others. Consequently, each Lagrangian model 
yields a distinctly different value of the GMR energy for a given value of K^. 
The model NLcr delivers the lowest GMR energy for a nucleus whereas the 
model SIGO produces the highest value amongst the three models for any 
given value of K^. 

The GMR energy is shown to exhibit a sensitivity to the Lagrangian model 
employed. This behaviour stems from the differences in the implicit density 
dependence of the meson masses in various models, thus altering the density 
dependence of the nuclear interaction in the RMF theory from one Lagrangian 
model to the other. The sensitivity to the density dependence opens up another 
dimension in the landscape of the GMR energy and incompressibility. 

The multiplicity of the GMR energy for a given in the RMF theory renders 
the 'microscopic' method of extraction of using interpolation amongst 
forces as inapplicable. For this approach to work, however, it is necessary to 
'calibrate' the density dependence of the meson masses and correspondingly 
the nuclear interaction in the RMF theory. 

Using the liquid-drop type expansion of the finite nucleus incompressibility 
Ka-i it is shown that for a given value of i^oo, each Lagrangian model yields 
a markedly different value of the surface incompressibility K sur f. The model 
NLcr delivers the largest value of K sur f, whereas the model SIGO provides 
the smallest K sur f amongst the three Lagrangian models considered. Conse- 
quently, a different value of the ratio K^fjK^ emerges in different models. 
Different K sur f values arising in various Lagrangian models are shown to be 
largely responsible for a multitude of GMR energies for a given Koo. Thus, the 
response of the nuclear surface to compression is found to be dependent on 
the Lagrangian model employed. Additionally, the ratio K^f/K^ with the 
Lagrangian models NLcr and SVI emerges to be significantly larger than the 
ratio of ~ — 1 that is usually obtained with the assumption of scaling in the 
breathing-mode GMR [12]. 

The GMR response in the non-relativistic Skyrme approach is simple in con- 
trast. The Skyrme approach commonly delivers a GMR energy which depends 
primarily on the value of [T3]. This difference in the behaviour of the 
GMR energies between the RMF theory and the Skyrme approach can be 
attributed to the fact that the form of the density dependence in the Skyrme 
density functional is well prescribed at the outset in marked contrast to the 
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various forms of density dependence which have emerged in the RMF theory. 
In a forthcoming work, we shall shed light on the differences amongst various 
Lagrangian models, which produce a varied response to the GMR energy in 
the RMF theory shown in this work. Further investigations will also be nec- 
essary in order to be able to discern and formulate the appropriate density 
dependence in the RMF theory. 
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